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■ The volume operator plays a pivotal role for the quantum dynamics of Loop Quantum Gravity (LQG) . 
I It is essential in order to construct Triad operators that enter the Hamiltonian constraint and which 
• become densely defined operators on the full Hilbert space even though in the classical theory the triad 

l/^ I becomes singular when classical GR breaks down. 

■ The expression for the volume and triad operators derives from the quantisation of the fundamental 
I electric flux operator of LQG by a complicated regularisation procedure. In fact, there are two inequivalent 

volume operators available in the literature and, moreover, both operators are unique only up to a finite, 
^ I multiplicative constant which should be viewed as a regularisation ambiguity. 

bi). Now on the one hand, classical volumes and triads can be expressed directly in terms of fluxes and 

this fact was used to construct the corresponding volume and triad operators. On the other hand, fluxes 
can be expressed in terms of triads and therefore one can also view the volume operator as fundamental 
and consider the flux operator as a derived operator. 
. In this paper we mathematically implement this second point of view and thus can examine whether 

the volume, triad and flux quantisations arc consistent with each other. The results of this consistency 
analysis are rather surprising. Among other findings we show: 1. The regularisation constant can be 
uniquely fixed. 2. One of the volume operators can be ruled out as inconsistent. 3. Factor ordering 
ambiguities in the definition of triad operators are immaterial for the classical limit of the derived flux 
operator. 

The results of this paper show that within full LQG triad operators are consistently quantized. 
In this paper we merely present ideas and results of the consistency check. In a companion paper we 
supply detailed proofs. 
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1 Introduction 



The major unresolved problem in Loop Quantum Gravity (LQG) (see J for books and J for reviews) 
is the satisfactory implementation of the quantum dynamics which in turn is governed by the Wheeler - 
DeWitt quantum constraint |S1I1], also called the Hamiltonian constraint. The volume operator lUllHI plays 
a pivotal role for the very definition of the Hamiltonian constraint because with its help one can quantize 
triad functions which enter the classical expression for the Hamiltonian constraint. In particular, one takes 
advantage of a Poisson bracket identity between the triads and the Poisson bracket among the Ashtekar 
connection and the classical scalar volume function. Namely, one uses the methods of canonical quantisation, 
that is, the axioms of quantum mechanics, according to which Poisson brackets between classical functions 
are turned into commutators between the corresponding operators divided by ih, at least to lowest order in 
h. Quite surprisingly, the Hamiltonian constraint operator and similarly also length operators are then 
densely defined on the full kinematical Hilbert space of LQG although the classical triad becomes singular 
in physically relevant situations such as black holes or the big bang. 

While playing such a distinguished role for the most important open problem of LQG, the volume operator 
and thus the derived triad operators have never been critically examined concerning their physical correctness 
and mathematical consistency. By the first we mean that it has never been shown within full LQG that 
the volume operator has the correct classical limit with respect to suitably chosen kinematical semiclassical 
states, for example those constructed in |S|^. 

By the second we mean the following: The fundamental kinematical algebra 21 on which LQG is based is 
the Holonomy - Flux algebra 10 and its representation theory together with background independence 
leads to a unique kinematical Hilbert space Now classically the volume and triad can be written 

as limits of functions of the flux. To implement them at the quantum level, one has to go through a 
complicated regularisation procedure and to take the limit. It is surprising that the resulting operators 
are densely defined at all and in fact have a discrete spectrum because they are highly non - polynomial 
expressions as functions of fluxes. This is the payoff for background independence, since in background 
dependent formulations, such as the standard Fock representations these operators are too singular. On 
the other hand, since there is little experience with non Fock representations, it is not at all clear whether 
the corresponding operators have anything to do with their classical counterpart. In fact, there are at least 
two ambiguities already at the level of the volume operator: First of all, there are in fact two unitarily 
inequivalent volume operators [SI El which come from two, a priori equally justified background independent 
regularisation techniques. We will denote them by Rovelli - Smolin (RS) and Ashtekar - Lewandowski (AL) 
volume respectively for the rest of this paper. Secondly, both volume operators are anyway only determined 
up to a multiplicative regularisation constant Creg which remains undetermined when taking the limit, 
quite similar to finite regularisation constants that appear in counterterms of standard renormaisation of 
ordinary QFT. The ambiguity is further enhanced by factor ordering ambiguities once we consider triad 
operators. These ambiguities are parameterized by a spin quantum number i = 1/2, 1,3/2, ... 
In this paper we will be able to remove all those ambiguities by the following consistency check: As we 
mentioned above, the volume and triad can be considered as functions of the fluxes. But the converse is 
also true: The fluxes can be written in terms of triads and thus the volume. Is it then true that there exists 
a regularisation constant for the volume operator and a factor ordering of the flux operator considered as 
a function of the triad operator or volume operator such that the corresponding alternative flux operator 
agrees (at least in the correspondence limit of large eigenvalues of the volume operator) with the fundamental 
flux operator, independent of the choice of £7 This better be possible as otherwise the inescapable conclusion 
would be that the volume operator is inconsistently quantised^. 

^This has been achieved, so far, only within a certain approximation |2] which essentiaUy consists in replacing SU{2) by 
Uilf. The necessary calculations in the Non Abelean case are much more complictaed due to two reasons: 1. The spectrum 
of the volume operator is not available in analytical form and 2. the semiclassical analysis is calculationally more difficult. 
However, work is now in progress in order to fill this gap. 

^In contrast, the triad operator follows from the volume operator by the axioms of quantum mechanics namely that Poisson 
brackets be replaced by commutators divided by ih and therefore it is not possible that the source of a possible problem is in 
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We will be able to precisely answer this question affirmatively. In more detail we will show: 

1. The RS volume operator is inconsistent with the flux operator, the AL volume operator is consistent. 

2. Creg = 1/48 Can be uniquely fixed, there is no other choice which is semiclassically acceptable. Re- 
markably, this is precisely the value that was obtained in 6; by a completely different argument. 

3. The choice of £ plays no role semiclassically, in fact it drops out of the final expression for the alternative 
operator altogether. Therefore the afore mentioned factor ordering ambiguity is absent as far as the 
flux operator is concerned. 

4. There is yet one more ambiguity in LQG which is already present classically: Classically it is possible 
to consider the electric field either as a two form or as a pseudo two form. The corresponding sign of 
the determinant of the triad is then either encoded in the electric field or in the conjugate connection. 
One can take either point of view without affecting the symplectic structure of the theory. We will be 
able to show that one must consider the electric field as a pseudo two form, otherwise the alternative 
fiux operator becomes the zero operator! 

5. As expected, the alternative and fundamental fiux operator agree for all values of the Immirzi param- 
eter ^Ij, hence it cannot be fixed by our analysis which is good because it has been fixed already by 
arguments coming from quantum black hole physics |15j . 

6. The calculations in this paper make extensive use of certain advances in technology concerning the 
matrix elements of the volume operator ^^l- Thus, our calculations provide an independent check of 
[T^ as well. 

7. The factor ordering of the alternative flux operator is unique if one insists on the principle of mini- 
mality^ . 

These results show that instead of taking holonomies and fluxes as fundamental operators one could instead 
use holonomies and volumes as fundamental operators. It also confirms that the method to quantise the 
triad developed in ,3^ is mathematically consistent. 

This paper is organised as follows: 

In section two we review the regularisation and definition of the fundamental fiux operator for the ben- 
efit of the reader and in order to make the comparison with the alternative quantisation easier. 
In section three we eplain the idea on which the construction of the usual fiux operator is based on. 
In section four we derive the classical expression for the alternative fiux operator. 

In section five and six we discuss our results for the alternative flux operator and compare them with the 
corresponding results of the usual fiux operator. 
Finally in section seven we summarise and conclude. 

All the technical details and tools that are needed to perform this consistency check are providid in our 
companion paper |19j . 



the quantisation of the triad operator. 

''By this we mean that given a function / with self adjoint quantisation / and a multiphcation operator g for which is 
defined everywhere on the Hilbert space, we may always consider instead /' — {gfg~^ + g~^ f'g) /2 as the operator corresponding 
to / if it has selfadjoint extensions as well. By minimalistic we mean the choice g — 1. This issue is always present even in 
ordinary quantum mechanics however it is usually not mentioned because one usually deals with polynomials and g ^ 1 would 
destroy polynomiality. In GR the expressions are generically non - polynomial from the outset and thus polynomiality is not 
available as a simplistic criterion. However, we may still insist on a minimal number of such kind of factor ordering ambiguities. 
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2 Review of the usual Flux Operator 



In LQG the classical electric flux Ek{S) through a surface S is the integral of the densitised triad E'^ over 
a two surface S 

E^{S) = jElnl (2.1) 

s 

where is the conormal vector with respect to the surface S. In order to define a coresponding flux 
operator in the quantum theory, we have to regularise the classical flux and then define the action of the 
operator on an arbitary spin network funtion (SNF) T^^^^ : G'^^'^^' — > C, where G is the corresponding 
gauge group, namely SU{2) in our case, as the action of the regularised expression, denoted by El{S) in 
the limit where the regularisation parameter e is removed 

EkiS)T^,rn,n ■■= ihYim {^^(5), T^,^,^} . (2.2) 

Here the limit is to be understood in the following way: The Poisson bracket on the right hand side of eqn 
(|2.2j) is calculated by viewing T^j,^^ as functions of smooth connections. After taking e — > one extends 
the result to functions of distributional connections and thus ends up with an operator defined on the Hilbert 
space of LQG. 
Classically, we have 

dT - 

{mS),T^J,^^A{heiA)h,Ei,))} = Y: {^r>(/^e)..}^7^. (2.3) 

Here {he)AB denotes the S'f/(2)-holonomy. The regularisation can be implemented by smearing the two 
surface S into the third dimension, shown in figure ^ so that we get an array of surfaces St- The surface 
associated with t = is our original surface S 

mS):=Y^ I dtEuiSt). (2.4) 



t 



t = +e 
t = 
t = -e 




St 

Figure 1: Smearing of the surface S into the third dimension. We obtain an array of surfaces St labelled by the 
parameter t with t E {— e, +e}. The original surface S is associated with t — 0. 

In order to derive the action of the flux operator on an arbitrary SNF, we would have to analyse the Poisson 
bracket among the flux and every possible SNF. Fortunately, it turnes out that each edge belonging to the 
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Figure 2: Edges of type up, down, in and out with respect to the surface S. 



associated graph 7 of can be classified as (i) up, (ii) down, (iii) in and (iv) out with respect to the 

surface 5. (See figure |2l for a graphical illustration.) An arbitrary T^^^^ contains then a particular amount 
of edeges of each type. Accordingly, if we have the knowledge of the Poisson bracket among the flux and any 
of these types of edges, we will be able to derive the Poisson bracket among E'^ and any arbitrary T^^^^. 
The calculation of the regularised Poisson bracket can be found e.g. in the second reference of 1 . After 
having removed the regulator we end up with the following action of the flux operator on a SNF 

Ek{S)T,,rn,n-^l, ^ , (2.5) 

e£E{'y) 



where is related to the Pauli matrices by := —iak- The sum is taken over all edges of the graph 7 
associated with ^. The function e(e, S) can take the values {—1,0, +1} depending on the type of edge 

that is considered. It is +1 for edges of type up, -1 one for down and for edges of type in or out. 
By introducing right invariant vector fields X^, defined by {X^f){h) := ^/(e*'^'-'/i) 
action of the flux operator as 



, we can rewrite the 

i=0 



Ek{S)T^J,r^^^ = -il ^ e{e,S)XtT^j^... (2.6) 

ee-E(7) 

Note, the right invariant vector fields fulfill the following commutator relations [Xg,X|] = —lerstX^. By 
means of introducing the self-adjoint right invariant vector field := — |Xg, we achieve commutator 
relations for that are similar to the one of the angular momentum operators in quantum mechanics 
[yj, y/] = ierstY^. Therefore, we also can describe the action of -Ea,.(S') by the action of the self-adjoint right 
invariant vector field on T ^ ^ ^ 

Ek(S)T ^ ^ ^ = --il y e(e, 5)K*^ T . - (2.7) 

ee-B(7) 



3 Idea of the Alternative Flux Operator 

Our starting point will be the Poisson bracket of the Ashtekar-connection Aa and the densitised triad 
given by 

{Ai{x),Etiy)]=5'ix,y)6t6^ (3.1) 

which we take as fundamental. In order to go from the ADM-formalism to the formulation in terms of 
Ashtekar varibales, one uses a canonical transformation. There exist two possibilities of choosing such a 
canonical transformation that both lead to the Poisson bracket above. These two possiblities are 



I Ai = Ti + J sgn{det{e)) Ki , E^ = ^ekrse^^'el 

II J^a = Ti + jKi , E^ = lekrse^^^el el sgn(det(e)) 
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Here, Ta is the SU(2)-spin connection, Kab = Kie'i, the extrinsic curvature (when the Gauss constraint 
holds) and 7 the Imirzi-parameter. Recah again the definition of the regularised classical flux E^{S) in eqn 
(|2.4|) . Now the idea of defining an alternative regularised flux 



2e 



dtEk{St); Ek{St 



Etn 



St 



(3.3) 



St 



is to express the densitised triad E'j: in terms of the triads as above. Due to the two possibile canonical 
transformations, we have also two possibilities in defining an alternative densitised triad 



rpa 



d et(e)eg 



abc„r 



2^k' 

1 , ^abc 

2 ^krs 
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b 



e""" sgn(det(e)) e\ e; 
S 



E\ 

e: 



k 

•a, II 



(3.4) 



where ea is the cotriad related to the intrinsic metric as q, 



ab 



eiel- From now on we will use E"^'^ and E[ 



a,ll 
k ' 



respectively for the two cases. 

So, instead of quantising the densitised triad directly, we could use the above classical identities, quantize 
them via the Poisson bracket identity in eqn ()4.1|) and check whether both quantisation procedures are 
consistent. 

The main difference between these two definitions is basically a signum factor which we will denote by S. 
From the mathematical point of view both definitions in eqn (|3.3|) are equally viable, thus we will keep both 
possibilities and emphasise the differences that occur when we choose one or the other definition. Notice 
however that dfii{E') = det(e)2 > gives an anholonomic constraint which appears to be inconsistent with 
the definition of Ek{S) as a derivative operator as in eqn 1)2. 5|) . as for instance pointed out in We will 
see that this is indeed the case. 

Inserting the alternative expression for into eqn 1)2. 4|) . we obtain 



EkiSt) = { 



abc^r „s„St 
b ^c"'a 



Er = det(e)e^ 



k 



(3.5) 



4 Construction of the Alternative Flux Operator 

The strategy for quantising the alternative Flux Operator will be as follows. As a first step we will apply the 
Poisson bracket identity in order to replace the triads eJJ by its associated Poisson bracket among the connec- 
tion Al and the scalar volume funtion V{R) = d^xy/det{q) = Jp, d-^x^/\ det(-E^)| = Jj^d^x^/\ det(£'")|. 
The Poisson bracket identity for the two cases is shown below 



{At,ViR)} 



K S 

2^b 



E\ 



k ~ 

<a,li 



y-a bc pS pt 



V_ abcCpS t 



(4.1) 



This relation is different for the two cases, because in deriving this relation we have to use the definition 
of the densitised triad in terms of the triads which is different for case I and case II. The difference 
between E'^'^ and E^^'" is again a signum factor S. Going back to eqn H3.5() we note that there is no S in 
E'^'' . Since we have to replace two triads e^,e^ by Poisson brackets, we get two factors of S in the case of 
S^'^. As is always one classically (since Qab is non-degenerate), it drops out in this case. In contrast, 
for E'^'^' we have a signum factor S occuring in the alternative flux in eqn ()3.5() . but no S in the Poisson 
bracket identity. Accordingly, we get only one S here, that does not dissapear classically, because it can 
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take the (constant) values ±1. Thus 



EiiSt) = I ekrse'''"'{AlViR)}{Al,Vmn^' 



Ei'iSt) = ^ I ekrse''''^{AlV{R)}S{Al,V{R)}n^,^ 



(4.2) 



St 



When later on we replace the classical expression by their corresponding operators, the main difference 
between Ej^{St) and E'f/{St) will be a so called signum operator S. Before, we have to replace the connections 
A^ by holonomies since A^ cannot be promoted to well defined operators in the Ashtekar-Lewandowski 
Hilbert space TIal- Hence, we choose a partition Vt of each surface St into small squares of area e'^. In the 
limit where e' is small enough we are allowed to replace the connection A^ along the edge by its associated 
holonomy h(ei). The partition is shown in figure ©• Usually this is done for holonomies in the fundamental 




Figure 3: Partition Vt of the surface St into small squares with an parameter edge length e'. 

representation of 1/2. But, as we want to keep our construction of the alternative flux operator as general 
as possible and to study the effect of factor ordering ambiguities we will consider holonomies with arbitrary 
representation weights i. The corresponding relation between the connection integrated along the edge e^, 
denoted by A^j from now on, and the associated holonomy is given by 



1 



{A^j{n), F(i?,c))} -TT.iTr) + o(e'2) = +vr,(/ie,) {Mh~'), V{R,ia))} 



(4.3) 



whereby we indicate a representation with weight £ by vr^ . 

Considering eqn (|4.3|) . we end up with the following classical identity for '•^^El{St) 



^'W.iSt) = hm efc--2 {AS(n),y(i?,c))} {^f(D),y(i?,c))} 

Devt 



lim > 



16 



1 



Tr 



(^"/(□))'^(«^C))})(4-4) 
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and in the case of '•'''El/{St) with 

= hm j;e,,,4{^3(°),mC))}5{vi:(n),y(i?.c))} 

nevt 

lim 



Tr (^TTtiKsin)) [-^e 7r^(rfc)57r^(/ie4(n)) ^Tri{h^^^^^), l^(i?„(n))|)i4.5) 

Here we used Tr (vr^(Tr.7r£(rfc))7r£(rs)) = — + 1)(2^ + 1) that is derived in appendix A of jT^ . 
Now, in LQG there exist two volume operators, one introduced by Rovelh and Smohn in 1994 (Vus) ,5, and 
another one pubhshed in 1995 by Ashtekar and Lewandowski (Val) 0- Hence, we have actually in each 
case two different possibilities for the Poisson bracket, because we could either use Vrs or Val- Thus, case 
I as well as case H splits into two different versions of alternative fluxes 

'''Ei'{St) - ''%'^'''{St), [St) (4.6) 

From now on we will use the notation above for the four different cases. Before we apply canonical quanti- 
sation, we want to discuss the two volume operators and their differences a bit more in detail. 

4.1 The two Volume Operators in LQG 

4.1.1 The Volume Operator Vrs of RoveUi and Smolin 

The idea that the volume operator acts only on vertices of a given graph was first mentioned in ^7]. The 
first version of a volume operator can be found in ^5, and is given by 



R 



= E J\^Crese,,kXl^Xi,X!^^^ . (4.7) 



I,J,K 

Here we sum over all triples of edges at the vertex v G V{'y) of a given graph 7. V^s is not sensitive to the 
orientation of the edges, thus also linearly dependent triples have to be considered in the sum. Moreover, 
we introduced a constant Creg £ that we will keep arbitrary for the moment and that is basically fixed 
by the particular regularisation scheme one chooses. When working with the volume operator we want to 
select physically relevant gauge invariant states properly. Hence, it is convenient to express our abstract 
angular momentum states in terms of the recoupling basis. The following identity jl2j holds 

^e,,,Xl^XlX^^ = \[Y,l,Y,\] =: \q^,,, (4.8) 

where Yu := Y^^ + Y^^ and Y^^ denotes the self-adjoint vector field Y^^ := -^X^^. 
Consequently, we get 



veV{'y)nRI<J<K 



V{R)^n JM;M'> = il Yl E V ^Creggfj^ \ JM; M'> . (4.9) 



yRS 
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The additional factor of 3! is due to the fact that we sum only over ordered triples I < J < K now. The 
way to calculate eigenstates and eigenvalues of V is as follows. Let us introduce the operator Ql]fjK 

Qv,IJK ■ ^p~7 

-CregqJjK (4.10) 

As a first step we have to calculate the eigenvalues and corresponding eigenstates for Q^\fjK- If for example 
\(t)> is an eigenstate of Qy'jjK with corresponding eigenvalue A, then we obtain V\(f>>= ^y\X\\4> >• Con- 
sequently, we see that while QZ]fjK can have positive and negative eigenvalues, V has only positive ones. 
Furthermore, if we consider the eigenvalues it A of Qv'jjk and the corresponding eigenstate |(/>+a>, |0_a>, 
we notice that these eigenvalues will be degenerate in the case of the operator V, as y^l + A| = y^l — A|. 

4.1.2 The Volume Operator Val of Ashtekar and Lewandow^ski 

Another version of the volume operator which differs by the chosen regularisation scheme was defined in [H] 



V{R)};^^'^\JM; M'> = £ >J J ^Creg >J e{e„ej,e^)qjj^ \ J M ; M'> . (4.11) 




The major difference between Val and Vrs is the factor e{ej,ej,eK) that is sensitive to the orientation of 
the tangent vectors of the edges {ej^ej^en}. e{ej,ej,eK) is +1 for right handed, —1 for left handed and 
for linearly dependent triples of edges. In the case of Val it is convinient to introduce an operator Q^"*^ 
that is defined as the expression that appears inside the absolute value under the square root in V^^ 

QV' := if^Creg ^ ' ^Lk (4-12) 

I<J<K 

By comparing eqn (|4.9|) with (|4.11j) we notice that another difference between V^s and Val is the fact that 
for the first one, we have to sum over the triples of edges outside the square root, while for the latter one, 
we sum inside the absolute value under the square root. Apart from the difference of the sign factor, the 
difference in the summation will play an important role later on. 



4.2 Factor Ordering Ambiguities 

The discussion in this section is valid for all four different version of the alternative flux, hence we will 
neglect the labels I,II,RS,AL here. 

The meaning of the limit when the regularisation parameter tends to zero, i.e. e — > in combination with 
the limit of the partition Vt when the edge length parameter e' — > is explaind in detail in in section 
4.3. It is the same limit as taken for the usual flux operator sketched between eqn ()2.2() and eqn 1)2. 3() . We 
have e.g. for those S which intersects 7 only in an edge of type up schematically 

Ek{S)T = hm ^ [ dtY,<T^\Ek{St)\T> T/ 

—e ^ 
''^ I 

where T/ are the SNW states that contribute to Ej^{ST)T. < T/ | Ek{St) | T > is actually independent of e 
and T/ — > Tg as a function of smooth connections which can then be extended to distributional ones. 
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It turns out that the alternative flux operator can be regularised to such an extent that each plaquette 
of the partition Vt of a surface St has an intersection with only one single edge of a graph 7 associated 
with an arbitrary given SNF T^^rn,n- If the considered edge lies completely inside the plaquette of St or 
completely outside (up to sets of dt measure zero if inside) it will lead to a trivial action of the alterna- 
tive flux operator. This is in full agreement with the usual flux operator. The alternative flux operator 

''^'£'fe(S't) attaches two additional edges, namely 63,64 to the edge e. The discussion in ^H] showed that 

these additional edges have to lie inside the surface St since otherwise the ^'^E^iSt) would have only again 

a trivial action, see also figure Furthermore, by attaching the additional edges 63, 64, cre- 



ei(n) 





flSt 






' e' 





64(0) /' \ 

/ 1 

62(0) 



Figure 4: An non-vaishing contribution to the matrix element of '^"'Ek on a given arbitrary SNF Tg, i.e. < 

Ts' I *'^'i?fe(n) I Ts > can only be achieved if Tg contains edges of type up and/or down, respectively with respect 
to the surface St- Moreover, the edges 63(0), e4(n) have to be attachted to Tg in this specific way. 

ates a new vertex at a certain point of e and thus divides e into an edge 61 of type up and an edge 62 
of type down with respect to St"^- Accordinly, as for the usual flux operator, the action of the alternative 
flux operator is totally determined by the action on edges of type up and down with respect to the surface St- 

Since we are familiar now with the action of '•^^EkiSt), we are able to consider the impact of different 
factor orderings. Going back to the classical identity in eqn H4.4() and eqn H4.5() respectively and considering 

the results of the discussion above, we know that ^'^E^iSt) adds two additional edges 63,64 to a given edge 
6 and devides this edge into an up edge 61 and an edge of type down 62 This is illustrated in figure El Let 
us call these SNF that involves the edges 61, 62 only l/J-'i^ ni2>. Here we use the so called recoupling basis to 
express the SNF and ji2 denotes the total angular momentum to which the two edges 61 and 62 couple at 
their single vertex v, whereas 71-12 is the associated magnetic quantum number. The SNF that results from 

1/5-^^^ ni2> by the action of ^'^'Ek^totiSt) and that contains four edges {61, 62, 63, 64} will be denoted by |a/ M>. 
In this case J is the total angular momentum, M the associated magentic quantum number and i as an 
index is needed, because with four edges more than one state exists with the same total angular momentum, 
but different intermediate couplings. For more details see section 5.2 and 6 in [E]- Now, dealing still with 
the classical expression in eqn ()4.4() and eqn (|4.5() respectively, we can apply the trace and rearrange the 
terms in a certain manner, because classically holonomies commute in order to obtain a sensible operator 
ordering in the quantum theory. If we consider ^^^Ej^'^'' (St) and ^^^E'f/'^^{St) that contain Val we know that 
due to the sign factor e{ej, 6.7, e^) in V^l the action of Val on linearly dependent triples vanishes. Therefore, 
it has to be ensured that the two holonomies ■ni{he^)EA and TTe{he^)iG act before Val does. This restriction 
reduces the number of possible factor orderings down to a single one. 

The situation is different for '*'£^^'^'^(5f) and ^'Wj^''^'^ (St), because here Vrs is involved, that has an non 

■^The opposite case is also possible of course, but we will restrict our discussion to the first case and emphasize in the following 
disscussion where exactly the choice of ei as a down edge and 62 of type up will make a difference. 
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Figure 5: The SNF \f3^^^ ni2> is transformed into an new SNF |af m> by the action of <"£'fc,t„t(S't). 



trivial action on linearly dependent triples. Consequently, more than one possible factor ordering exists. Let 
us discuss the number and differences of these orderings later on and restrict ourselves for ^'^Ej;^^ (St) as well 
as ^'^EjJ''^'^ (St) to the single ordering that is possible for Val, first. It turnes out that for all operators with dt 

measure 1, the dependence of Ek{St) on t drops out and taking the average ^ f dt becomes trivial resulting 

— e 

in eqn 1)4. 13() . We will thus keep the label St in what follows, but keep in mind that the t-dependence is 
trivial. 



4.3 Canonical Quantisation 

Usually the densitised triads, appearing in the classical flux Ek{S) are quantised as differential operators, 
while holonomies are quantised as multiplication operators. If we choose the alternative expression -Efc(S') we 
will instead get the scalar volume V and the so called signum S operator into our quantised expression. The 
properties of this S will be explained in more detail below. Moreover, we have to replace Poisson brackets 
by commutators, following the replacement rule {. , .} — > {l/ih)[., .]. The detailed derivation of the final 
operator can be found in section 4 of citeGT. 

Clearly, we want the total operator to be self-adjoint, so we will calculate the adjoint of '^'^Ek{St) and define 

the total and final operator to be ^'^Ek,tot{St) = ^{Ek{St) + -E'^(S't)) that is self-adjoint by construction. 
Hence, the final operator for V^s which we will use through the calculation of this paper is given by 



^I/II.RS 

'''E,.o. (St) 



-'^e{Tk)cB'^e{(^)Ei 



ez)FG 



whereby we used the identity Tti{hf,^)AB = ^t{hei)BA]\ T^e{^)AB = (—1)^ "^^a+b^o and the definition of the 
Planck length := {hK)~'^. The box surounding the signum operator S should indicate that it is contained 
in the operator in case II, while it is not in case I. 

Considering the operator Val, we know that for each commutator only one term will contribute, because 
otherwise we cannot construct linearly independet triples of edges since {ei, 62, 63/4} are linearly dependent. 
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Therefore in the case of Val we obtain the following final expression 



'■e^ )eA 



Here again for case II the signum operator is inlcuded, whereas in case I it is not. 
Next, we want to calculate the matrix elements < /^■'^^ rn-ia |*^'i?fc,tot(5't)| rriia > of all four versions 

^I,AL ^I,RS ^11, AL ^U.RS 

'''E,^,JSt), ''^E,^,JSt), (St) and ^'^E,^,^, (5*) of the new flux operator. 

The action of the holonomy operators ■ni{hej)AB on \(3^'^'^ mi2> can be described in the framework of angular 
momentum recoupling theory with the powerful tool of Clebsch-Gordan-Coefficients (CGC). The correspon- 
dence between the Ashtekar-Lewandowski Hilbert space TCal and the abstract angular momentum Hilbert 
space is discussed in section 5.1 in JHl- Hence, the matrix element will roughly speaking have the following 
structure 

oc C{j,2j]d3ma-^)C%j,2j]a.,ma^)C%d.,J-,JM)C{a-,,i-,JM) < M\d\a-j M> (4.16) 

Here C{ji,j2] J M) denotes the CGC that we get if we couple the angular momenta ji and j2 to a resulting 
angular momentum J with corresponding magnetic quantum number M and O is an operator including 
the operator Vr^ and Val respectively and 5. The details can be found in section 5.1 of jT^]. Due to the 
symmetry properties of the alternative flux operator that are analysed in section 5.2 of |19j . the only total 
angular momenta that will contribute to the final matrix element are J = 0, 1. (See section 5.2 of |19j for 
more explanations concerning this point.) Since the physically relevant states are gauge-invariant we choose 

ji2 = 0. The behaviour of ^^^Ek^tot{St) under gauge transformations, that is discussed in section 5.3 of ^Hl 
leads to the restriction jij = 1- Consequently, at the end of the day we get the following form of the matrix 

element of '■'^^Ek^totiSt) 



<(3^ m,2\^''>Ek,USt)\(i^o> 

8£-4(_i)3^ +^ . 

i2+B+F,0 

< 1 ; £ -B I ^ mi2+B> < £ fhrz+B ; £ F I 0> 

< a^M = m-ia+B+F ; m' \ 0-i\a\M = ; m' m'^y 

-(—1) '^F+B,0<Jc+F,mi2 

< 00;£C|£C> < £C;If\IC+F> 

+ ^ < OL^M = ■,'fn[ m'2 \O2\01\M = c+F ; m' 

\/2lTT 1 ^ ^, ^, , ^ , 1 , , 

j=z — < a^M = mi2 \rn^ | C2 M = C+F ; rn^> 



V3 



< a\M = mj^2\'rn[m'^\02\a\M = c+F ■,m[m',^> |, (4.17) 
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whereby the four different cases are encoded in the operators Oi and O2. Explicitly, we have 



^1 — ^ AL 

^1 ~ ^134 ~'~ ^234 ~'~ ^134^234 + ^234^134 + ^134^123 + ^124^134 + ^234^123 + ^124^234 + ^124^123 

^2 ~ ^ AL 

^2 ~ + ^234 ~'~ ^234^134 + ^134^234 + ^123^134 + ^134^124 + ^123^234 + ^234^124 + ^123^124 

Of '^^ = ValSV^l 

of '-^^ = +^9134^^5134 + ^g234'^V'q234 + ^"9134^^9234 + ^9234^^9134 + ^9134^^9123 + ^9124^^9134 + ^9234*^^9123 
+ ^124*^^234 + ^124*^^123 



Of '^■^ — +^9134^^9134 + ^9234'^V'q,234 + ^9234^^5134 + ^"9134^^9234 + ^9123^^9134 + ^9134^^9124 + ^9123^^9234 

+ + ^91235^^9124, (4.18) 

whereby we used the notation Vq^jj^ for Vrs = Vq^^^ + Vq^^^ + Vq^.^^ + Vq^2'i when only the triple {ej, Cj, e^} 
contributes to Vrs. The derivation of the various versions of Oi,02 can be found in ^Hl- These 8 versions 
result from taking into account a) the two volume operators, b) the signum operator 5 or not and c) the 
adjoint or not in eqn (|4.17|) . 



Before we will discuss our results in the next sections, we want to explain a bit more in detail what we 

mean by say the consistency check is affirmative or not. We managed to implement an alternative flux 

^ ^i,al ^i,rs ^ii.al ^ii,rs 

operator ^'^k^USt) in four different versions ''^Ej,^,^,{St), ^'^E,^,JSt), ^'^E,^,^, (St) and ^'^E,^,^, (St). The 
usual flux operator is quantised as a differential operator in the standard way. The alternative flux operator 

^^^Ek,tot{St) is quantised via the Poisson bracket identity in eqn (|4.1|) analogous to the quantisation of the 
Hamiltonian Constraint. If these two methods of quantisation are mathematically consistent with each 

other, the action of Ek{S) and the one of *^'-E'fc,tot('S'i) should only differ by a constant, namely 



Ek{S)\p\>= C{j,e)Cre/'E,^,JS)\p\>= C{j,£)CregY. < '''E,^,JS) |/?%> 1/3^ mi2>, (4.19) 

mi2 

^0 ^ 

whereby '•'^Ej^ ^^^{S) denotes the alternative flux operator ^'^^Ek,tot{S) where Creg has been replaced by 1. 
The constant C{j,i) might depend on the spin labels j,l of the edges, but at least semiclassically the 

dependence on the spin label j and £ must disappear since otherwise the behaviour of Ek{S) and ^'^Ek,tot{S) 
in the correspondence limit of large j would disagree. Notice that also a dependence of lim C{j,i) =: C{£) 

on i is inacceptable, because classically the flux is independent of the factor ordering ambiguity i. Moreover 
this constant will fix the ambiguity in the volume operator Creg that is due to regularisation as we will see 
later. Hence, lim C{j,i) = 1/Creg = const as will be discussed more in detail in section IHH 



^I,RS -^I,AL 

5 Case I: Results for ^^^E.^^^^St) and ^^E.^^^^St) 

^I,AL 

5.1 Calculations for ^"^fc,,„,(5t) 

For technical reasons, we consider only a spin label i = 0.5, 1, because higher spin labels cannot computed 
analytically anymore. Fortunately, the main properties of this case already occurs when considering small 
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L The detailed calculation in section 6.2 in ^H] show 

< = ■,fh[ m'2 1 0['^'" \aiM = c+F ; m[ m'^> = 

< aj M = ■,fh[ \ Ol;"^^ \a\M = c+F ; m' m'^y = 0, (5.1) 

where i = 2,3,4. Going back to eqn (|4.17|1 . we note that the vanishing of the matrix elements above has 
the consequence that the whole matrix element < /J-^^^ m-ia I'^'-E^ tot{St)\ Z?-'^^ 'm^2> is zero. Since the action 
of A; toi(<S't) on an arbitrary SNF can be derived from exactly this matrix element, we can conlcude that 
''^^^k,tot{^t) is the zero operator. Accordingly, ^^^Ef,^^^^{St) is not consistent with the usual flux operator. 

5.2 Calculations for <''^fc ^„,(5't) 

Also for the operator '^'^Ef, t^X^t) the analogous calculations discussed in section 6.3 of ^Hl yields only trivial 
matrix elements 

< = fh'^ \ 0['^^ |a? M = C+F ; m[ m'^y = 

< a\ M = ifh'-^ rh'2 1 O^'^^ \ a\ M = c+F ; m' m'2> = 0, (5.2) 

Consequently, we can draw the same conlcusion as for E ^^^{St) and state that '•^^ E f. ^^^{St) is the zero 

operator and therefore inconsistent with the usual flux operator. Futhermore as either E ^^^[St) nor 

'^'-^fe toti^t) survive the consitency check, we can rule out, at least for the cases of ^ = 0.5, 1, the choice of 
E^ = det(e)e^ on which these operators are based on. To rule out the choice E'^{St) = det(e)e^ completely, 
we need to investigate the matrix element for arbitrary representation weights i. For higher values of i the 
calculation cannot be done analytically any more. However, the results for i = 0.5, 1 indicate that there is 
an abstract reason which leads to the vanishing of the matrix elements for any i. We were not able to find 
such an abstract argument yet. However, even if that was not the case and there would be a range of values 
for I for which not all of the matrix elements would vanish, it is inacceptable that the classical theory is 
independent of i while the quantum theory strongly depends on £ in the correspondence limit of large j. 

6 Case II: Results for ^^^E,^^^^ (St) and ^^^E,^^^^ {St) 

^II,AL 

6.1 Calculations for {St) 

<:ii,AL ^ 

Considering the case of the operator '•^-'Ej.^^^ (St), we can read of from eqn (|4.1<S|) the expressions Oi = 
ValSVal = 02- Since the signum operator S that corresponds to the classical expression S := sgn(det(e)) 
does not exist in the literature so far, we will in detail explain how the operator S has to be understood. 

6.2 The Signum Operator S 

We are dealing now with case II meaning that the densitised triad is given by E'^'" = 5det(e)e^, where 

ia,'- 
k 

det(S) = sgn(det(e))det(g) with det(g) = [det(e)]^ > 0. (6.1) 



5 := det(e). Applying the determinant onto E'^'" , we get 



Therefore, we obtain 

sgn(det(^)) = sgn(det(e)) = S. (6.2) 
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In the following we want to show that S = sgn(det(i?)) can be identified with the signum of the expression 
inside the absolute value under the square roots in the definition of the AL-volume. For this purpose let us 
first discuss this issue on the classical level and afterwards go back into the quantum theory and see how 
the corresponding operator S is connected with the operator Q^^ in eqn ()4.12|) . 

In order to do this let us consider eqn (|4.5j) . This equation contains the classical volume where 

denotes a region centred around the vertex v{d). 
The volume of such a cube is given by 



^y\det{E)\cfx, 



(6.3) 



R 



R. 



where we used det(g') = | det(i?)| from eqn 1)6.11) . Introducing a parametrisation of the cube now, we end up 
with 







I 


duj 







y/\det{E){u)\d^u= J \det{X)\ ^/\d^t{E)(uj\d^u. 



(6.4) 



2 2 J 



In order to be able to carry out the integral we choose the cube Ry(ci) small enough and thus, the volume 
can be approximated by 



,'3 



det(-)(.) 



V\det{E){v)\. 



(6.5) 



Using the definition of det(E') = ^Cabc^-^^^'EjE^Ef, we can rewrite eqn (|6.3|) as 



V{Kiu)) = J ^l^eabce^'^^E-ElEf 



(6.6) 



If we again choose -R^(n) small enough and define the square surfaces of the cube as , we can re-express 
the volume integral over the densitised triads in terms of their corresponding electric fluxes through the 
surfaces S' 



1 

3!' 



The flux through a particular surfaces S' is defined as 

E,{S') = I E'^nf nf -- 



e^kiE,{S')Ek{SJ)Ei{S'<) 



(6.7) 



(6.5 



Here nf^ denotes the conormal vector associated with the surface S' . Regarding eqn (|6.7j) we realise that 
inside the absolute value in eqn 1)6. 7() appears exactly the definition of det(£^j(S'^)). Therefore we get 



(6.9) 



On the other hand, by taking advantage of the fact that the surfaces S' are small enough so that the integral 
can be approximated by the value at the vertex times the size of the surface itself, we obtain for det(£'j(5^)) 

dQi{Ej{S')) ^ det{E^{v)n^\v)e'^) 



det(£;|(z;))det(nf (z;))e' 
det{E{v))det{nf{v))e'^ 



(6.10) 
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If we consider the definition of the normal vector in eqn (|(i.8j) ,we can show the following identity 

det(nf ) = det{Xf det{X-^) = = det{Xf. (6.11) 

Inserting eqn (|6.11|) back into eqn (|6.1U|) we have 

detiEjiS')) « det{E{v))[det{X{v))]h'^ (6.12) 

and can conclude that eqn (|6.9j) is consistent with the usual definition of the volume in eqn (|6.5|1 . 
Since we want to identify S := sgn(det(i?)) with the signum that appears inside the absolute value under 
the square root in the definition of the volume, we can read off from eqn llOl) . that we still have to show 
sgn(det(i?)) = sgn(det(£'j (5^))). However, this can be done by means of eqn H6.12|) 

sgn{det{Ej{S'))) ^ sgn{det{E{v))[det{X{v))fe'^) 

= sgn{det{E{v))) sgn([det(X(w))]2) sgn{e'^) 

= sgn(det(^(i;))). (6.13) 

Consequently, we must identify S with the signum that appears inside the absolute value under the square 
root in the definition of the volume V in the classical theory, because it was precisely the expression 
det{Ej{Si)) that was used in the construction of the volume operator, defined as the square root of absolute 
value of det{E). In the quantum theory, we introduced the operator Q^^ in eqn (|4.12|) . which is basically 
the expression inside the absolute value in the definition of the volume operator. Hence, it can be seen as 
the squared version of the volume operator that additionally contains information about the signum of the 
expression inside the absolute values. Consequently, we have the operator identitity Q^^ = ValSVal- Now 
we will be left with the task to calculate particular matrix elements for Q^^ which can be done by means 
of the formula derived in ■ 

6.3 Calculations for ^''^^ [St) 

One big advantage that comes along with the operator identity = ValSVal is that diagonalisation of 
the operator is no longer necessary as it was in case I for V^^, because the operator '•'E^, {St) contains 
only particular matrix elements of Q^^ that can be exactly calculated, even for arbitrary I, by means of 
the tools developed in The details of this calculation can be found in |19j as well as the corresponding 
matrix elements of the usual fiux operator Ej.{S). If we compare the results of the usual flux operator with 

^11, AL ^11, AL 

the one of '•^^E;. (St) we can judge whether ^'^Ej. (St) leads to a result consistent with the usual fiux 
operator. It transpires 

^n,AL ^ 

(5)1/3° o>= 3!8C,e,i^fc(5)|/30 o> (6.14) 

Therefore the two operators differ only by a positive integer constant. As there is still the regularisation 
constant Creg in the above equation we can now fix it by requiring that both operators do exactly agree 
with each other. In fact there is no other choice than exact agreement because the difference would be a 
global constant which does not decrease as we take the corresponding limit of large quantum numbers j. 
Thus, we can remove the regularisation ambiguity of the volume operator in this way and choose Creg to be 

r — J- - J- 

<^reg ■— 3I8 — 48- 

This is exactly the value of Creg that was obtained in [H] by a completely different argument. Thus the 

geometrical interpretation of the value we have to choose for Creg is perfectly provided^. 

Note that the consistency check holds in the full theory and not only in the semiclassical sector. Conse- 

^11, AL 

quently, the operator '•'^Ej. (St) is consistent with the usual flux operator. 

^The factor 8 = 2^ comes from the fact that during the regularisation one integrates a product of 3 5— distributions on R 
over R"'' only. The factor 6 = 3! is due to the fact that one should sum over ordered triples of edges only. 
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^11, RS 

6.4 Calculations for (St) 

^11, RS 

Now, considering the operator '•^''Ej^ (St) things look differently. Here, a quantisation that is consistent 
with Vjis of the signum operator S cannot be found due to the simple reason that Vrs in contrast to Val is 
a sum of single square roots. Hence, there is no origin for a global sign as it was in the case of Val- (See 

^I,RS 

also section 6.6.1 in JH].) In retrospect there is a simple argument why the only possibility '•'^ E ^^^{St) 

^11, RS 

(since ^^'-E^ {St) does not exist) is ruled out without further calculation: Namely, the lack of a factor of 

orientation in Vjis, like e(ej,ej,ex) in Val, leads to the following basic disagreement with the usual flux 
operator: Suppose we had chosen the orientation of the surface S in the opposite way. Then the type of the 
edge e switches between up and down and similarly for 61,62. Then, the result of the usual flux operator 
would differ by a minus sign. In the case of Val we would get this minus sign as well due to 6(67,6^,6^) 
contained in Q^' = ValSVal, whereas a change of the orientation of 6^, 62 would not modify the result of the 
alternative flux operator if we used Vrs instead, because it is not sensitive to the orientation of the edges. 

Accordingly, we should stop here and draw the conclusion that '^'-E^ i^t) is not consistent with Ei.{S). 

^ ^ ^ ^11. RS ' ^ 

One might propose to artificially use V^sSalVrs for ^'Ef. (St). Note, that we attached the label AL to S to 
emphasize that its quantisation is in agreement with Val- This is artificial for the following reason. Suppose 
we have a classical quantity A := det{E) and two different functions /i := y^|^| and /2 := sgn(^). If we 
want to quantise the functions /i and /2, we do this with the help of the corresponding selfadjoint operator 

A and obtain due to the spectral theorem /i = ^J \A\ and /2 = sgn(A). The product of operators VhsSalVus 
rather corresponds to gi = A' and g2 = sgn(^), because V^s is quantised with a different regularisation 

scheme than S is. This would only be justified if ^|J\A\ and A' would agree semiclassically. However they do 
not: If we compare the expressions for Val and V^s then, schematically, they are related in the following way 
when restricted to a vertex: = I x^'-es E eiei, ej, eK)qijK\^^'^ while Vy^ns = E l^CregQuKl^^'^ - 

I<J<K I<J<K 

It is clear that apart from the sign e(e/,6j,6^) the two operators can agree at most on states where only 



one of the quK is non vanishing (three or four valent graphs) simply because yj\a + b\ 7^ + vH for 
generic real numbers o, h. 

r:ii,RS ^ ^ ^ 

Nevertheless by analysing ''^'-E^ {St) when the artifical operator VhsSalVrs is involved, we obtain 

^11. RS ^ 
'^^E;,,,,, (S)|/3%>= C{j,i)CregEkiS)\(3\>, (6.15) 

whereby C{j,i) G R is a constant depending non-trivially on the spin labels j,i in general. One can show 
that C{j, £) C{tj semiclassically, i.e. in the limit of large j, which is shown in appendix E and discussed in 

^II.RS ^ ^ ^ 

section 6.6.2 of of JHI- Hence ''^^E^. (St), including the artificial operator Vrs^alVrs, would be consistent 
with Ei.{S) within the semiclassical regime of the theory if we chose Creg = ^/C{£) and if C{£) would be 
a universal constant. Unfortunately, C{i) has a non-trivial i -dependence which is inacceptable because it 
is absent in the classical theory. Moreover, we do not see any geometrical interpretation available for the 
chosen value of Creg for any value of i in this case. One could possibly get rid of the ^-dependence by simply 
cancelling the linearly dependent triples by hand from the definition of V^s- But then the so modified Vj^g 
and Val would practically become identical on 3— and 4— valent vertices and moreover V^g now depends on 
the differentiable structure of E. See more about this in the conclusion section. 

7 Conclusion 

We hope to have demonstrated in this paper that at least certain aspects of LQG are remarkably tightly 
defined : Certain factor ordering ambiguities turn out to be immaterial, some regularisation schemes can be 
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ruled out as unphysical once and for all. The fact that we can exclude the RS volume from now on as far 
as the quantum dynamics is concerned should not be viewed as a criticism of jSj at all: The regularisation 
performed in is manifestly background independent, natural and intuitively very reasonable. It was the 
first pioneering paper on quantisation of kinematical geometrical operators in LQG and had a deep impact 
on all papers that followed it. Most of the beautiful ideas spelled out in 'Sj also entered the regularisation 
performed in ^ and continue to be valid. One could never have guessed that the regularisation performed 
in [2] leads to an inconsistent result. It took a decade to develop the necessary technology in order to 
perform the consistency check provided in this paper. Hence, the fact that we can define the theory more 
uniquely now should be taken as a strength of the theory and not as a weakness of iSj- Interestingly, the 
very first paper on the volume operator |12| that we are aware of does look more like Val rather than Vns- 
We speculate that if one had completed the ideas of one would have ended up with Val rather than 
Vns- Of course, one could take the viewpoint that the consistency check performed here is unnecessary, that 
one can just take some definition of the volume operator and not worry about triads. However, as triads 
prominently enter the dynamics of LQG such a point of view would render the quantum dynamics obsolete. 
In other words, the dynamics and all other operators which depend on triads such as the length operator 
or spatially difeomorphism invariant operators forces us to use Val rather than Vas- Finally notice that 
one of the motivations for choosing V^s rather than Val is that V^s does not depend on the differentiable 
structure of S. Hence one could use homeomorphism rather than diffeomorphism in order to define "diffeo- 
morphism invariant" states as advertised in |21j.|2Uj. This makes the Hilbert space of such states separable. 
However, notice that homeomorphisms are not a symmetry of the classical theory. Furthermore, there are 
other possibilities to arrive at a separable Hilbert space: One can decompose the Hilbert space Ti-Diff corre- 
sponding to the strict diffeomorphisms into an uncountably infinite direct sum of separable Hilbert spaces. 
In the current proposals for the quantum dynamics all of these mutually isomorphic Hilbert spaces are left 
invariant. If the theory are left invariant by the Dirac observables, then they would be superselected and 
any one of them would capture the full physics of LQG. 

The technical details of the check, mainly displayed in our companion paper 19^, required a substantial 
computational effort. There are an order of ten crucial stages in the calculation where things could have 
gone badly wrong. Following the details of our analysis, one sees that all the subtle issues mentioned must 
be properly taken care of in order to get an even qualitatively correct result. These subtleties involve, among 
other things: 

1. The meaning of the limit as we remove the regulator and to define the alternative flux operator has 
to be understood in the same way as for the fundamental flux operator, otherwise the alternative flux 
opertor is identical to zero. This issue is discussed in full detail in section 4.3 of |19j 

2. Switching from the spin network basis to the abstract angular momentum basis, discussed elaborately 
in section 5.1 of JH]) has to take care of the precise unitary map between these two representations of 
the angular momentum algebra, otherwise the two operators differ drastically from each other. This 
unitary map is not mentioned in the literature because for gauge invariant operators it drops out of 
the equations. However, for the non - gauge invariant flux it has a large impact. 

3. Very unexpectedly, the sign of the determinant of the triad enters the calculation in a crucial way: 
Classically one would expect it to be neglegible, especially in an orientable manifold. However, had 
we dropped it from the quantum computation then the alternative flux operator would again vanish 
identically. It is very pleasing to see that quantisation based on a pseudo vector density is ruled out. 
This is for the same reason that one cannot implement the momentum operator ih-^ on ^2(11''', dx). 

4. A different ordering than the one we chose would have resulted again in the zero operator even in case 
II. 

5. As for the fundamental flux operator one has to first smear the alternative flux operator into the 
direction transversal to the surface under consideration. For the fundamental flux this implies that 
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edges of type "in" are not acted on. Without this additional smearing the fundamental flux would 
be ill defined. For the alternative flux this implies furthermore that the classification of edges into 
the types up, down, in and out is meaningfull at all. Indeed, one actually can define the alternative 
flux without the additional smearing. The result is well defined. However, it would differ drastically 
from the fundamental flux as soon as there is a vertex of valence higher than two of the graph of the 
spin network state in question within the surface. The additional smearing has the effect that with dt 
measure one all the vertices within the sufaces St are bivalent. 

6. Furthermore, the analyses in section 4.3 of 'W shows that without the additional smearing we would 
be missing the crucial factor 1/2 in eqn (|4.13|) and our Creg would be off the value found in [£]. 

7. Just following the tedious calculations step by step as explicitly shown in jj^l and evaluating all the 
CGC's involved one sees that all the subtle signs have to be there, everything fits only when doing the 
calculation with 100% accuracy. The calculation is therefore a highly sensitive consistency check. 

8. All the i dependence disappears even at small values of j. This is especially surprising because the 
classical approximation of a connection by a holonomy along a given path becomes worse as we 
let £ grow. Of course in the limit we take paths of infinitesimal length, however, this is done after 
quantization and it could have happened that the quantization is affected by a non trivial i dependence 
which however should disappear in the limit of large j. 

It transpires that the reason for getting a zero operator without the signum operator unveils a so far not 
appreciated symmetry of the volume operator. It would be desirable to understand the symmetry from a 
more abstract perspective. 

This paper along with our companion paper is one of the first papers that tightens the mathematical 
structure of full LQG by using the kind of consistency argument that we used here. Many more such checks 
should be performed in the future to remove ambiguities of LQG and to make the theory more rigid, in 
particular those connected with the quantum dynamics. 

Acknowledgements 

It is our pleasure to thank Johannes Brunnemann for countless discussions about the volume operator. 
We also would like to thank Carlo Rovelli and, especially, Lee Smolin for illuminating discussions. 
K.G. thanks the Heinrich-Boll-Stiftung for financial support. This research project was supported in part 
by a grant from NSERC of Canada to the Perimeter Institute for Theoretical Physics. 

References 

[1] C. Rovelli, "Quantum Gravity", Cambridge University Press, Cambridge, 2004 

T. Thiemann, "Modern Canonical Quantum General Relativity", Cambridge University Press, 2005; 
'ir-qc/0110034 

[2] C. Rovelh, "Loop Quantum Gravity", Living Rev. Rel. 1 (1998) 1, gr-qc/9710008 

T. Thiemann, "Lectures on Loop Quantum Gravity", Lecture Notes in Physics, 631 (2003) 41 - 135, 
'i^-qc/0210094 

A. Ashtekar, J. Lewandowski, "Background Independent Quantum Gravity: A Status Report", Class. 
Quant. Grav. 21 (2004) R53; |i7^c/0404018 

L. Smolin, "An Invitation to Loop Quantum Gravity", hep-th/0408048 
[3] T. Thiemann, "Anomaly-free Formulation of non-perturbative, four-dimensional Lorentzian Quantum 
Gravity", Physics Letters B380 (1996) 257-264, gr-qc/9606088 

T. Thiemann, "Quantum Spin Dynamics (QSD)", Class. Quantum Grav. 15 (1998) 839-73, 
gr-qc/9606089 ; "II. The Kernel of the Wheeler-DeWitt Constraint Operator", Class. Quantum Grav. 



19 



15 (1998) 875-905, "gr-qc/9606090'; "III. Quantum Constraint Algebra and Physical Scalar Product 
in Quantum General Relativity", Class. Quantum Grav. 15 (1998) 1207-1247, gr-qc/9705017 ; "IV. 
2-|-l Euclidean Quantum Gravity as a model to test 3+1 Lorentzian Quantum Gravity", Class. Quan- 
tum Grav. 15 (1998) 1249-1280, |gr-qc/9705018| ; "V. Quantum Gravity as the Natural Regulator of 
the Hamiltonian Constraint of Matter Quantum Field Theories", Class. Quantum Grav. 15 (1998) 
1281-1314, gr-qc/9705019 

[4] T. Thiemann, "The Phoenix Project: Master Constraint Programme for Loop Quantum Gravity", 
|gr-qc/0305080 

B. Dittrich, T. Thiemann, "Testing the Master Constraint Programme for Loop Quantum Gravity I. 
General Framework", gr-qc/0411138 ; "II. Finite Dimensional Systems", gr-qc/0411139 ; "III. SL(2,R) 
Models ", |gr-qc/041114 "IV. Free Field Theories", gr-qc/0411141 ; "V. Interacting Field Theories", 
|gr-qc/041 11421 



[5] C. Rovelli, L. Smolin, "Discreteness of volume and area in quantum gravity", Nucl. Phys. B442 (1995) 

593, Erratum : Nucl. Phys. B456 (1995) 734 
[6] A. Ashtekar, J. Lewandowski, "Quantum Theory of Geometry II : Volume Operators", Adv. Theo. 

Math. Phys. 1 (1997) 388-429 
[7] T. Thiemann, "A Length Operator for Canonical Quantum Gravity", J. Math. Phys. 39 (1998) 3372- 

3392; gr-qc/9606092 

[8] T. Thiemann, "Quantum Spin Dynamics (QSD): VII. Symplectic Structures and Continuum Lat- 
tice Formulations of Gauge Field Theories", Class. Quant. Grav.l8:3293-3338, 2001, hep-th/0005232"; 
"Gauge Field Theory Coherent States (CCS): I. General Properties", Class. Quant. Grav. 18:2025- 
2064,2001, hep-th/0005233 

T. Thiemann, O. Winkler, "Gauge Field Theory Coherent States (GCS): II. Peakedness Properties", 
Class.Quant.Grav.l8:2561-2636,2001, hep-th/0005237 ; "III. Ehrenfest Theorems", Class. Quantum 
Grav. 18 (2001) 4629-4681, hep-th/0005234 ; "IV. Infinite Tensor Product and Thermodynamic Limit", 
Class. Quantum Grav. 18 (2001) 4997-5033, hep-th/0005235 

H. Sahlmann, T. Thiemann, O. Winkler, "Coherent States for Canonical Quantum General Relativity 

and the Infinite Tensor Product Extension", Nucl.Phys.B606:401-440,2001; gr-qc/0102038 

T. Thiemann, "Complexifier Coherent States for Canonical Quantum General Relativity", 



|gr-qc/02 6037 



[9] H. Sahlmann, T. Thiemann, "Towards the QFT on Curved Spacetime Limit of QGR. 1. A General 
Scheme", [gr-qc/0207030 ; "2. A Concrete Implementation", "gr-qc/0207031" 
[10] H. Sahlmann, "When do Measures on the Space of Connections Support the Triad Operators of Loop 
Quantum Gravity?", gr-qc/0207112 ; "Some Comments on the Representation Theory of the Algebra 
Underlying Loop Quantum Gravity" , |gr-qc/0207111| 

H. Sahlmann, T. Thiemann, "On the Superselection Theory of the Weyl Algebra for Diffeomorphism In- 
variant Quantum Gauge Theories" , ,gr-qc/0302090 ; "Irreducibility of the Ashtekar- Isham-Lewandowski 
Representation" , _gr-qc /0303074| 

A. Okolow, J. Lewandowski, "Diffeomorphism Covariant Representations of the Holonomy Flux Alge- 
bra" , gr-qc/0302059" 

C. Fleischhack, "Representations of the Weyl Algebra in Quantum Geometry" , "math-ph/0407006| 
J. Lewandowski, A. Okolow, H. Sahlmann, T. Thiemann, "Uniqueness of Diffeomorphism Invariant 
States on Holonomy - Flux Algebras" , to appear 
[11] A. Ashtekar, C.J. Isham, "Representations of the Holonomy Algebras of Gravity and Non-Abelean 
Gauge Theories", Class. Quantum Grav. 9 (1992) 1433, hep-th/9202053 

A. Ashtekar, J. Lewandowski, "Representation theory of analytic Holonomy C* algebras", in "Knots 
and Quantum Gravity", J. Baez (ed.), Oxford University Press, Oxford 1994 



[12] T. Thiemann, "Closed Formula for the Matrix Elements of the Volume Operator in Canonical Quantum 

Gravity", J. Math. Phys. 39 (1998) 3347-3371; gr-qc/9606091 
[13] R. Loll, "Imposing det E ^ in discrete quantum gravity," Phys. Lett. B 399 (1997) 227 

|arXiv:gr-q c/9703033'. 




20 



[14] G. Immirzi, "Quantum Gravity and Regge Calculus", Nucl. Phys. Proc. Suppl. 57 (1997) 65; 
gr-qc/970 1052. 

C. Rovelli, T. Thiemann, "The Immirzi Parameter in Quantum General Relativity", Phys. Rev. D57 
(1998) 1009-14; [gr-qc/9705059 
[15] A. Ashtekar, J. C. Baez, K. Krasnov, "Quantum Geometry of Isolated Horizons and Black Hole En- 
tropy", Adv.Theor.Math.Phys.4:l-94,2001, gr-qc/0005126 

K. Meissner , "Black Hole Entropy in Loop Quantum Gravity", Class. Quant. Grav. 21 (2004) 5245- 
5252; [i^qc/0407052 

M. Domagala, J. Lewandowski, "Black Hole Entropy from Quantum Geometry", Class. Quant. Grav. 21 

(2004) 5233-5244, 2004; g r-qc/040 7051 
[16] J. Brunnemann, T. Thiemann, "Simplihcation of the Spectral Analysis of the Volume Operator in Loop 

Quantum Gravity", gr-qc/0405060 

[17] L. Smolin, "Recent Developments in Non-Perturbative Quantum Gravity", hep-th / 9202022 



[18] J. Brunnemann, private communication 
[19] K. Giesel, T. Thiemann, "Consistency Check on Volume and Triad Operator Quantisation in Loop 

Quantum Gravity 11", gr-qc/0507037 
[20] W. Fairbairn and C. Rovelh, J. Math. Phys. 45 (2004) 2 802 arXi v:gr-qc/0403047_. 
[21] N. Grot and C. Rovelh, J. Math. Phys. 37, 3014 (1996) |arXiv:gr-q c/9604010 . 



21 



